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ABSTRACT 


Various transition methods are used here to study the viscous effects encoun- 
tered in low density, hypersonic flight, through the transition from free molecular 
to continuum flow. Methods utilizing Viking Data, Shuttle orbiter data, a Potter 
number parameter and a Shock Reynolds number were implemented in the Program to 
Optimize Simulated Trajectories (POST). Simulations of the Aeroassist Flight 
Experiment (AFE) using open loop guidance were used to assess the aerodynamic 
performance of the vehicle. A bank angle was found for each transition method that 
would result in a 200 nautical-mile apogee. 

Once this was done, the open loop guidance was replaced by the proposed guidance 
algorithm for the AFE. Simulations were again conducted using that guidance and the 
different transitions for comparison. For the gains used, the guidance system showed 
some sensitivity in apogee altitude to the transition method assumed, but the 
guidance was able to successfully complete the mission. 
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INTRODUCTION 

Hie development of an orbiting Space Station has generated a need for a space - 
based, reusable vehicle, capable of transferring large payloads from a high energy 
orbit to a Low Earth Orbit (LEO). Although this type of maneuver can be done propul- 
sively through a Hohmann transfer, the amount of propellant required limits the pay- 
load capacity. Another method is an aeroassisted maneuver vhich utilizes the aerody- 
namics of the vehicle by dipping into the Earth's upper atmosphere and expending 
energy through drag and heat, ‘ Once the vehicle exits, a comparatively smaller 
change in velocity is needed to circularize the orbit at apogee, thus providing 
significant propellant savings. Remaining in the atmosphere too long, however, will 
result in extreme heat rates, and excessive energy loss vhich will prevent the 
vehicle from escaping. 

Because the success of the mission depends heavily on the energy loss induced by 
the atmospheric pass, the aerodynamic performance of the vehicle is important. From 
previous missions it can be seen that in low density, hypersonic flight, vehicles 
experience a large decrease in performance due to viscous effects encountered in the 
transition from free molecular to continuum flow. 1 J Rarefied, hypersonic flow 
causes: (1) a large increase of friction effects; (2) a moderate increase of some 

pressures; (3) wall and shock slip; and (4) a thickening and merging of the shock and 
boundary layers. Current research facilities are unable to simulate the conditions 
encountered in this type of situation, and the existing data are insufficient for 
further AOTV design studies. Hie Aeroassist Flight Experiment (AFE) is an effort to 
obtain further data that will aid in the verification of the computational codes 
needed to design the AOTV's. 

The AFE will be deployed by the Shuttle orbiter in LEO and will fire its Solid- 
propellant Rocket Motors (SRM) in order to give it the velocity necessary to simulate 
an AOTV return from geosynchronous orbit. After passing through the atmosphere and 
expending enough energy to achieve the required target apogee, the vehicle will con- 
tinue to apogee where it will circularize to rendezvous with the Shuttle (Figures 1 
and 2). In order to predict realistically the trajectory of the mission for instru- 
ment calibration purposes and guidance development, the viscous effects encountered 
in the transition regime must be taken into account. Because no agreement exists 
on a fundamental parameter to correlate accurately the effects on aerodynamic 
performance, only the continuum values have been used previously. Four existing 
methods are studied and discussed in this paper: (1) Viking data correlation; ^ 

(2) Lockheed bridging formula; 4 (3) a Potter number; and (4) the Shock Reynolds 

number transition. [6] Because none of the methods have been verified in a direct 
flight application, the guidance algorithm used in the guidance system of the AFE 
must be universal enough to be able to handle a reasonable amount of error. 

Modifications were made in the three degree-of-freedom version of the Program to 
Optimize Simulated Trajectories (POST) to include tw subroutines: (1) one with four 

viscous calculation options (Appendix A); and (2) one with the current proposed 
guidance algorithm for the AFE. In order to rate the aerodynamic performance of the 
AFE, open loop guidance versions of the program were run with the four different 
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aerodynamic transition models, and with the standard continuum values to be used for 
comparison. Closed loop guidance versions using the proposed algorithm were also run 
in order to test its ability to handle the different bridging models and still 
reasonably meet the targeted conditions.^ Hie same initial conditions and event 
criteria were used for all the runs and are shown on Figures 1 and 2. 

TRANSITION FROM FREE MOLECULAR TO CONTINUUM FLOW 

The first opportunity to study the effects of low density hypersonic flight on a 
vehicle entering the atmosphere arose with the entry of the Mars probe, the Viking 
lander, into the Martian atmosphere. Using flight measurements from pressure 
instruments, accelerometers, and a mass spectrometer, the previously unknown values 
of the drag coefficient between the free molecular and continuum values were 
defined. C3] Using these data, Jim Jones of the Langley Research Center established 
an equation which calculates values for the drag coefficient in the transition regime 
as a function of a parameter called VBAR. VBAR is the Mach number divided by the 
square root of the free stream Reynolds . number. 

Another method was derived using the Shuttle Aerodynamic Design Data Book. The 
Lockheed Bridging Formula that comes from these data requires only two endpoints, 
which are the free molecular and continuum values for the drag coefficient. This 
formula is modeled after the relatively high pressure drag component of hot cylin- 
ders. The rarefaction effects of the transition hypersonic flow are accounted for 
using a fundamental parameter called the Knudsen number. This parameter is equal to 
the free stream mean free path divided by the characteristic length of the vehicle, 
which is the diameter in this case. This bridging formula is unique to the design of 
the vehicle, and is based solely on empirical data. In the studies done in this 
paper the coefficient of lift is also varied using the Lockheed Bridging Formula. 

The third method to be discussed in this paper was formulated by Leith Potter of 
Vanderbilt University, who is an expert in low density fluid dynamics. The parameter 
used here. Potter number, is a Reynolds number corrected for the enthalpy of the 
flow. An argument proposed by Potter is that lifting flight through different types 
of flow involves changes in angles of attack and even changes in the shape of the 
vehicle. To account for these effects, he includes an SSTAR in his Potter number 
calculation which is the cross sectional area divided by the watted area of the 
vehicle raised to the one half power. Sometimes there is a problem with calculating 
the wetted area when the lines of separation are hard to determine. A variety of 
drag coefficient data obtained from spheres, blunt nosed cones, lifting bodies and 
the Shuttle orbiter closely correlate with the Potter number. [5] 

The fourth transition method involves a correlating parameter called the Shock 
Reynolds number (R^). Ihe Shock Reynolds number is defined as 

p V D 
^00 00 


where is the free stream density, V^, is the free stream velocity, D is a 

reference length, and is the viscosity behind the shock. [6] p or r S 2 < i and 
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Mach number >14, lift, drag, and moment coefficients tend to correlate with R e ^ 
for blunt bodies. A representation of C D and C L variations with R 02 between 
free molecular flow and continuum flow has been developed and these transition equa- 
tions have been included in this investigation. 

All four transition methods resemble an exponential curve fit between the free 
molecular and continuum hypersonic flight drag coefficients. Figure 3 shows the dif- 
ference between the two drag coefficient values versus angle of attack, and Figure 4 
compares the transition methods used to bridge the gap. Figure 5 shows the gap be- 
tween the two lift coefficient values versus angle of attack, and Figure 6 compares 
the coefficient of lift variation used in the Lockheed Bridging and R 02 methods 
with the constant value used in the Viking and Potter versions. It should be. noted 
that the aerodynamics of the vehicle were cut off once the vehicle reached 400,000 
feet, which is considered to be the edge of the atmosphere. 


OPEN LOOP AERODYNAMIC SIMULATIONS 

Before a guidance algorithm could be implemented into the program it was neces- 
sary to determine the vehicle's performance capabilities. Comparative trajectories 
using the four available transition lift and drag coefficients, along with the con- 
tinuum lift and drag coefficient were examined using POST. A target apogee of 200 
nautical miles, and a heat rate not exceeding 180 Btu/ft^-sec was required for the 
mission to be considered successful. It should be noted that the heating rates shown 
here do not include corrections for viscous effects like those included for the 
aerodynamic coefficients. Therefore, this is not a valid evaluation of the actual 
heating rates. The vehicle was kept at a constant angle of attack of 1 7 degrees, and 
the program continued to alter the bank angles until the mission objectives were 
met. The bank angle vfoich worked for the standard version was then used in the other 
three versions in order to determine the sensitivity of the maneuver. The results of 
these runs can be studied in Table I. It can be seen that a change of less than 
0.1 degrees in the bank angle can result in an error of almost 8 nautical miles in 
the projected apogee. This can also be interpreted as a very small change required 
to put the AFE back on course in case something unforeseen arises. Overall the AFE 
has the aerodynamic capability necessary to make this mission a success. What is 
needed is a guidance algorithm vhich can use these performance capabilities to their 
best advantage. 


CLOSED LOOP GUIDANCE SIMULATIONS 

The guidance algorithm implemented in POST was developed to meet the demands for 
an aerobraking trajectory guidance technique that was uncomplicated, easily inte- 
grated and adaptable to a range of vehicle aerodynamic configurations. Included in 
the requirements was the ability to handle dispersions in entry conditions, atmo- 
spheric conditions, and aerodynamic characteristics. This algorithm uses the bank 
angle to control the lift vector so the vehicle will retain only the magnitude of en- 
ergy needed at the exit point in order to achieve the target apogee. 

Roll reversals are used to control the inclination at the exit point. No 
orbital plane change is wanted in, this case, so the vehicle must bank left and 
right. The roll rate is limited to a maximum of 15 degrees per second in all the 
cases presented here, and the reaction control system is not activated so the angle 
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of attack stays constant. Other than the the differences described above the same 
conditions were used here as in the open loop runs. 

During the initial use of the guidance routine, some unsteady behavior in the 
bank angle was observed. Discussion with the authors of reference 5 indicated that 
this unsteady behavior had been reduced by changing the method of maintaining the 
desired inclination angle. This new guidance algorithm was used in this study. ^ 

The guidance version was run for all four transition cases, and again it must be 
noted that no viscous corrections were included for the heat rate calculation. The 
results are shown in tabular form in Table II and on plots in Figures 7-16. In all 
four cases the mission objectives were achieved with an error of less than .01 
degrees in the inclination, but the apogee attained varied by as much as 8 nautical 
miles. Errors of up to 20 nautical miles in apogee can be compensated for so these 
errors were not excessive. Hie change in velocity required at apogee to circularize 
the orbit does not exceed the budget allocation. Ihe heat rate, minimum altitude, 
and maximum acceleration all fell within the acceptable design limits of the vehicle. 
Ihe bank angles commanded by the guidance are shown as Figures 17 through 21 for the 
five transition methods. As can be seen the character of the command is affected by 
the transition assumed. 


CONCLUSIONS 

Plots of the transition values for the drag coefficient were similar and seem to 
bridge the gap reasonably. 

The C D using the Potter number transition never reaches the continuum value, 
emphasizing the effects included in the SSTAR parameter. 

The trajectory parameters did not vary greatly regardless of vhich transition 
method was used in the simulation and the guidance algorithm essentially handled the 
mission in each case. However, the transition method assumed changed the character 
of the commanded bank angle and the apogee altitude showed some sensitivity to the 
transition method. 
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APPENDIX A 


"SUBROUTINE VISAERO" 
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Table I. Transition Method Comparison with Aerodynamic Model 


Transition 
Method 
Used for 

C L and C D 

$ for a 
200 NM 
Alta 

AV to 

Circularize 
at Apogee 

Max Heat 
Rate 

Btu/f t^-sec 

Max 
Accel 
(g's ) 

Minimum 

Altitude 

Inclination 

Angle 

400,000 ft. 

Projected 
Alta Using 
<j>=97. 367 (NM) 

Continuum 
C L and C D 

97.367 

334.079 

155.794 

2.598 

249 065 

32.246 

200.000 

Viking 
Data C D 

97.420 

334.733 

1 55.758 

2.598 

245 423 

32.242 

191 .527 

Lockheed 

Bridging 

C D and C L 

97.416 

334.536 

1 55.785 

2.568 

245 419 

32.238 

192.090 

Potter 
Number C D 

98.427 

337.418 

154.347 

2.609 

245 749 

32.125 

85.753 

Shock 

Reynolds 

Number 

C L and C D 

98.1 07 

338.300 

1 54. 734 

2.573 

245 628 

32.1 61 

86.430 


Table II. Transition Method Comparison with Guidance Model 


Transition 
Method 
Used for 
C L and C D 

AV to 

Circularize 
at Apogee 

Max Heat 
Rate 

Btu/f t^-sec 

Max 
Accel 
(g's ) 

Minimum 

Altitude 

Inclination 

Angle 

400,000 ft. 

Apogee 

Achieved 

Continuum 
C L and c d 

322.80 

148.448 

2.274 

1 

28.50 

198.0 

Viking 
Data Cq 

318.20 

148.378 

2.272 

251 829 

28.50 

195.0 

Lockheed 

Bridging 

c o and C L 

320.10 

148.407 

2.273 

248 010 

28.50 

195.9 

Potter 
Number C Q 

320.4 

146.433 

2.264 

248 670 

28.49 

192.1 

Shock 

Reynolds 

Number 

C L and C D 

315.4 

154.8 

2.58 

245 583 

28.50 

189.2 
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Initial Flight 
Characteristics of Proposed 
AFE 


Weight 

Velocity 

Sref 

Gamma 

Alpha 

L 

ref 


2741 lbs. 
3382.6 ft/sec. 
154 sq. ft. 
-4.5 deg. 

17 deg. (trim) 
14 ft. 


1976 Standard Atmosphere assumed 


Figure 1. AFE basic structure - propulsion module, and heat shield with a 
thermal protection system. 
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AFE Mission Events 

1 . Enter Atmosphere at 
400000 ft. 

2. Exit Atmosphere at 
400000 ft. 

3. Circularize at 200 Nm 
Apogee . 


Figure 2. AFE trajectory profile. 
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Figure 4 


Comparison of' drag coefficient from several transition methods versus 
time and at 17° angle of attack. 




Figure 5. Free molecular and continuum lift coefficient uersue angle of attack 
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Figure 6. Lockheed, Shock Reynolds, and continuum lift coefficient versus time. 
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Figure 7. Time histories of inertial velocity, altitude, and heat rate for runs 

with standard C and C and guidance active. 
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CD TIME, SEC 

CL 

Figure 8. Time histories of drag coefficient, lift coefficient, and altitude 


for rune with standard lift coefficient and drag coefficient and guidance 
active. 
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Figure 9. Time histories of velocity, altitude, and heat rate for runs with 
Viking lift coefficient and drag coefficient and guidance active. 
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Figure 10. Time histories of drag coefficient, lift coefficient, and altitude 
for runs with Viking drag coefficient and guidance active. 
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Figure 11. Time histories of velocity, altitude, and heat rate for runs with 

Lockheed lift coefficient and drag coefficient and guidance active. 
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Figure 12. Time histories of drag coefficient, lift coefficient, and altitude 
for runs with Lockheed drag coefficient and lift coefficient and 
guidance active. 
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Time histories of velocity, altitude, and heat rate for runs with 
Potter lift coefficient and drag coefficient and guidance active. 
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Figure 14. Time histories of drag coefficient, lift coefficient, and altitude 
for runs with Potter drag coefficient and guidance active. 


23 




HEATRT, 

BTU 

2 

FT SEC 

1130 


50 | 

I 

I 

! 

I 

i 

0 — 


i 

0 . a — | 

ALTITO, 

NM. 

0.6— i 

i 

I 


0 . 4 — | 

I 

j 

0 . 2 — 1 


UEL I 

RLTITG 

HERTRT 



TIME, SEC 


Figure 15. Time histories of velocity, altitude, and heat rate for runs with 

Shock Reynolds lift coefficient and drag coefficient and guidance active. 
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Figure 16. Time histories of drag coefficient, lift coefficient, and altitude 
for runs with Shock Reynolds lift coefficient and drag coefficient 
guidance active. 
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Figure 17. Bank angle time history for no transition. 
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Figure 19. Bank angle time history for Lockheed transition. 
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Figure 21. Bank angle time history for Shock Reynolds transition. 
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